Vortex Dipole in a BEC with dipole-dipole interaction 
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We consider single and multiply charged quantized vortex dipole in an oblate dipolar Bose Einstein 
condensate in the Thomas-Fermi (TF) regime. We calculate the critical velocity for the formation of 
a pair of vortices of opposite charge. We find that dipolar interactions decrease the critical velocity 
for a vortex dipole nucleation. 
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C ! I. INTRODUCTION 



A vortex dipole is a pair of vortices of equal and opposite circulation situated symmetrically about the origin. Under 
linear motion of a localized repulsive Gaussian potential, a vortex pair formation with opposite circulation is possible 
—I : if the potential is moved at a velocity above a critical value [ij . In the experiments [2j 01 j ^ laser beam focused 
^ on the center of the cloud was scanned back and forth along the axial dimension of the cigar shaped condensate. 
^jTj. Vortices were not observed directly, but the strong heating only above a critical velocity was measured. It was shown 
I ' that the measurement of significantly enhanced heating is due to energy transfer via vortex formation Recently, 
(-H . experimental observations of singly and multiply charged quantized vortex dipoles in a highly oblate BEC with 

■ were reported by Neely et al. (^J. In the experiment, vortex dipoles were created by forcing superfluid around a 
^ ] repulsive Gaussian obstacle using a focused blue-detuned laser beam. The beam was initially located on the left of 
C7^, the trap center and the harmonic potential was translated at a constant velocity until the obstacle ends up on the right 
of the trap center. At the same time, the height of the obstacle is linearly ramped to zero, leading to the generation of 
a vortex dipole that is unaffected by the presence of an obstacle or by heating due to moving the obstacle through the 
C ' edges of the BEC where the local speed of sound is small. Vortex dipoles were observed to survive for many seconds 
in the condensate without self-annihilation. The experiment also provided evidence for the formation of multiply 
(~| charged vortex dipoles. The authors in [5,] noted that the theoretical predictions of critical velocity for vortex pair 
O ] formation in Q are in good agreement with the experimental results. The critical velocity is given by the minimum 
O of the ratio of the energy to the momentum of the vortex dipole 

Vc = mm^^—) (1) 

\ where E{I) is the energy of an elementary excitation with linear momentum (or impulse) / [7[. When the object 
l/^ ■ moves at a velocity above a critical value, the superfluid flow becomes unstable against the formation of quantized 
' vortices, which give rise to a new dissipative regime Pairs of vortices with opposite circulation are generated 

at opposite sides of the object. Recently, instead of removing the trapping potential and expanding the condensate 
to make the vortex cores optically resolvable, Freilich et al. experimentally observed the real-time dynamics of vortex 
dipoles by repeatedly imaging the vortex cores [ll|. Vortex tripoles have also been observed exp erimentally Jl^. 
Several theoretical investigations have been reported for the generation ^13-17], stability and stationary 

configurations of vortex dipoles (20l . [2l| . In addition, fully analytic expressions of the angular momentum and energy 
of a vortex dipole in a trapped two dimensional BEC were obtained [22l | . 

The successful realization of Bose-Einstein condensation of ^^Cr atoms has stimulated a growing interest in the study 
of BEC with nonlocal dipolc-dipole interactions [23l - l29j . A particular interest is the vortex structures in dipolar 
condensates (30l - l46j . In contrast to the isotropic character of contact interaction, long ranged and anisotropic dipole- 
dipole interaction has remarkable consequences for the physics of rotating dipolar gases in TF limit. It was shown 
that, in axially symmetric traps with the axis along the dipole orientation, the critical angular velocity, above which a 
vortex is energetically favorable, is decreased due to the dipolar interaction in oblate traps [l^]. It was discussed that 
the effect of the dipole-dipole interaction is the lowered precession velocity of an off-center straight vortex line in an 
oblate trap [S] . Our aim in the present work is to calculate the critical velocity for vortex dipole formation in a dipolar 
oblate BEC in TF regime. This paper is structured as follows. Section II reviews Bose-Einstein condensates with 
dipole-dipole interaction. Section III investigates the critical velocity in the presence of the dipole-dipole interaction. 
The last section discusses the results. 
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II. DIPOLAR BEC 



Consider a BEC of N particles with magnetic dipole moment oriented in the z direction. In the mean field theory, 
the order parameter V'(r) of the condensate is the solution of the Gross-Pitaevskii equation (GPE) psi. [29j 

--^V2 + Vt+ 5l^(r)P + V-W = fii:{r) (2) 

2m / 

where /i is the chemical potential, g = as is the s-wave scattering length, Vr is the trap potential 

m 

Vt = \mujl (p2 + ^2^2) (3) 
7 is the trap aspect ratio, and $(j!i(r) is the dipolar interaction 

^ l\ 3 /" >3 /I - ^COS^fl 2 ,A\ 

•^ddiv) = —gsdd j d\' ^^_^,^3 W{v)\' (4) 

where r — r' is the distance between the dipoles, 9 is the angle between the direction of the dipole moment and the 
vector connecting the particles, and the dimensionless quantity Sdd is the relative strength of the dipolar and s-wave 
interactions 47[. The BEC is stable as long as —0.5 < Sdd < 1 in the TF limit l30l. l47l . In the regime Sdd < 0, the 



dipolar interaction is reversed by rapid rotation of the field aligning the dipoles [31|, |32 
We shall use scaled harmonic oscillator units (h.o.u.) for simplicity. In this system, the units of length, time, and 
I h 1 

energy are \ , — and hLu±, respectively. Hence the CP equation in h.o.u reads 

V 2m(jj± 2(jjj_ 

(-V^ + V^ + .g'|^'(r)p + ^'(r) = ^,'^'{r) (5) 
where the dimensionless interaction parameter g' is given by 



, 2mN 2muj^_ 

and ^'dd — ^ddig^g')- The normalization of i^'{r) chosen here is J d^r\^'{r)\'^ = 1. 

The equation ([S]) is an integro-differential equation since it has both integrals and derivatives of an unknown wave 
function. The solution of this equation was presented by Eberlein et al. in TF regime [47|. They showed that a 
parabolic density remains an exact solution for an harmonically trapped vortex-free dipolar condensate in TF limit: 

n(r)=no(l-^-^) (7) 

where hq = — and R and L are the radial and axial sizes of the condensate, respectively. The condensate aspect 
9' 

R 

ratio K is defined as k = — . In the absence of dipolar interaction, the condensate aspect ratio k and the trap aspect 

ij 



ratio 7 match, k decreases with increasing Sdd in an oblate trap. Hence, for Edd > [sdd < 0), k < 7 (n > 7) [30ll3l|. 
In h.o.u, the dipolar mean-field potential inside the condensate in the TF approximation is given by [431 

$,,(r) = nog'sddij^- ^ (l - 2W~L^)) 



where /(k) for oblate case, > 1), is given by 



2/ A arctan — 1 
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III. VORTEX DIPOLE 



Vortices can be nucleated in a BEC by a localized potential moving at a velocity above a critical value [48j . Vortices 
with opposite circulation are generated at opposite sides of the condensate. The motion of each vortex arises from its 
neighbor. In the presence of weak dissipation, the two vortices slowly drift together and annihilate when the vortex 
separation is comparable with the vortex core size. 

The critical velocity expression for a vortex dipole formation derived by Crescimanno et al. @ was found in good 
agrement with the experimental observation Q . We now extend their approach to include dipolar interaction. Con- 
sider a pair of single vortices with opposite charge. We suppose the vortices are located symmetrically about the trap 
center. Let the distance between the cores be represented by d. Substituting ?/'(r) — \J rt(r)e*"^ into the GP equation 
and equating imaginary and real terms leads to the following hydrodynamics equations: 

(^-V^ + (V0)2 + P1±^ + + Cl>,,^ = (10) 



- V^^ - 2V0-Vy^ = (11) 

The term V'^\fn in the equation (jlOp is negligible within TF approximation. The ansatz for the phase function for a 
vortex dipole is Q 



/ sin — ^ sin t , ^ 



0(r) = / ( arctan( ^^^ ) - arctan( '""" '^^ ) I (12) 

\ cos a cos B 



where p and 9 are the polar coordinates, and / is vorticity. The condensate velocity is given by v(r) — — V(/>(r). 

TO 

We assume that the vortices are far enough from each other, but near enough to the trap center. The ansatz is 
equivalent to requiring — < S < [i [61. Using (|12p we find 

|V0| = ( ^— I (13) 

' ' \^p2rf2cos2 0+(p2_ i|)2+^y ^ > 

where rj = . Excluding the vortex core regions from the domain complicates the analytic evaluation of energy 

and impulse precisely where the TF approximation fails. To prevent this difficulty, 77 is added to the denominator 
of (|13p [w. This regulated expression is confirmed by the observation that for vortex pair not too far from the trap 
center {d < fj,'), the contribution to y/n{r) from the kinetic energy term |V(/)p is never larger than /i' [6]. 
It is reasonable to approximate the TF density in h.o.u. by 

The correction due to the term |V(/ip is at the order of d^ /R*, which is very small compared to in h.o.u. Let us 
calculate the total energy and the total impulse of the condensate. The expression of total energy of the condensate 
in h.o.u is given by 

E = Jd^r (^(Vv^)2 + (^^V0)2 + (^!i_l!f!)^ + ^,,2 ^ ^^5) 

ih 

The total expression of impulse using the momentum of the condensate P ~ -^[(V'0*)V' ~ '0*Vi/'] is 



d^r|P| j d^'v n{v) |V0| (16) 
The critical velocity for vortex pair creation using the Landau criterion is defined as @ 



El — Eq 



(17) 
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FIG. 1: The critical velocity, Vc (mm/s), for Sdd = (solid) and Edd ~ 0.15 (dashed) as a function of d in units of h.o.u. for an 
oblate trap with 7 = 5 



— 1=1 

---1=2 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 



FIG. 2: The critical velocity, Vc (mm/s), for d = 1 in units of h.o.u. as a function of Edd for an oblate trap with 7 = 5. The 
solid (dashed) curve is for singly (doubly) quantized vortex dipole. 



Here Ei and // are the energy and impulse of vortex state whereas is the energy of non- vortex state. Analytical 
evaluation of these energy, impulse and critical velocity functions were performed by Crescimanno et al. in two 
dimensions for a nondipolar condensate [g*] . In this study, we will calculate the critical velocity of a dipolar condensate 
in an oblate trap. As the resulting equations for Vc are complicated, it is necessary to obtain them numerically for 
given as, Edd, We compare the critical velocities of dipolar and non-dipolar condensates. 



IV. RESULTS 



In this paper, within the TF regime we perform a numerical calculation of critical velocity for a vortex pair 
formation in a dipolar BEC. We examine dipolar gas containing 150000 ^^Cr atoms in an oblate trap with trap 
frequencies uj± = 2tt x 200 rad/s, uj^ =2tt x 1000 rad/s, so the trap aspect ratio is 7 = 5. The magnitude of the 
scattering length for ^^Cr is lOSas (a^ is the Bohr magneton). We assume that vortex separation d satisfies the 

1 2 / ■ • • ■ • 

condition — < d < /i . The chemical potential is approximately 42 in h.o.u. for Sdd ~ and changes very slightly 

. . . . 

with d. The dipole-dipole interaction decreases the chemical potential. Below, we perform numerical integration of 
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FIG. 3; The ratio of critical velocity to the speed of sound, — , for d = 1 in units of h.o.u. as a function of Edd for an oblate 

c 

trap with 7 = 5. The solid (dashed) curve is for singly (doubly) quantized vortex dipole. 
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energy and impulse functional (jl5ll6l) to find the critical velocity of vortex pair formation (jl7l) . 

It is well known that in an oblate trap, condensate aspect ratio, k, decreases with increasing Edd and the dipole-dipole 
interaction energy is positive. Minimizing the energy functional of the dipolar condensate, we find that the condensate 
aspect ratio increases slightly with d and the radius R is almost the same for all d. For example, k is between 4.73 and 
4.76 when d is between 0.5 and 6, respectively for Edd — 0.15. The kinetic energy changes significantly with vortices 
separation d, while the other terms in the energy expression change slightly with d. The impulse of vortex state, /, 
depends strongly on d and Sdd since it scales as d/R^ (|13ll6p . It decreases with Edd, since the dipolar interaction 
stretches the cloud radially in an oblate trap. 

The critical velocity Vc decreases with increasing separation d for a non-dipolar condensate @. This is because the 
TF density is a maximum at the trap center and reduces with the distance away from the trap center. We expect 
that the critical velocity decreases with d also for a dipolar BEC since parabolic form of density retains in the case of 
dipolar interaction. Fig-[T] plots the critical velocity Vc as a function of the distance between the vortices for Sdd = 
(solid curve) and for Edd — 0.15 (dashed curve). The critical velocity is between 1.25 — 2.15 mm/s for Edd — 0.15, 
and 1.27 — 2.24 mm/s for Edd = in the range 6 > d > 0.5. The critical velocity difference between dipolar and 
non-dipolar condensates becomes smaller as d is increased. This is because k increases with d. As can also be seen 
from the figure, the inclusion of dipolar interaction decreases the critical velocity for a fixed value of d. This is always 
true for positive values of Edd- In the case of negative values of Edd, the effect of dipole-dipole interaction is the 
increased critical velocity. Fig-[2] shows the critical velocity as a function of Edd for fixed d = 1. The solid curve shows 
singly quantized vortex dipole while the dashed curve shows doubly quantized vortex dipole. The effect that Vc is 
decreased with increasing Edd in an oblate trap is the first main result of this paper. It is energetically less expensive to 
nucleate a vortex in an oblate dipolar Bose-Einstein condensate than in a condensate with only contact interactions. 
At first sight, this might seem counterintuitive since the dipole-dipole interaction energy is positive in an oblate trap. 
It is remarkable to note that although dipolar interaction is positive for an oblate trap, the excess dipolar energy is 
negative. 

The nucleation of multiply charged vortex dipoles was observed for trap translation velocities well above Vc in the 
experiment 5]. Furthermore, it was observed that the vortices exhibit periodic orbital motion and vortex dipoles may 
exhibit lifetimes of many seconds, much longer than a single orbital period 5]. Fig-[2] compares the critical velocity 
for singly and doubly quantized vortex dipole. As expected, Vc is bigger for doubly quantized vortices. 
On the investigation of a vortex dipole, not only Vc, but also the ratio Vc/c is of importance. Here c is the speed of 
sound [131 



where a is the angle between directions of the wave vector and the dipoles. Suppose the direction of the phonon 
wave vector is perpendicular to the orientation of the dipoles (a = 7r/2). In this case, the speed of sound becomes 



c = w — g{l — edd)- Remarkably, both Vc and c decrease with Edd- However, the ratio Vc/c increases with Edd- As Edd 



goes to one, the critical velocity approaches to the speed of sound. In Fig-[31 we plot Vc/c versus Edd for singly (solid 
curve) and doubly quantized vortex dipoles (dashed curves) for fixed d — 1. The effect that Vc/c is increased with 
increasing Edd in an oblate trap is the second main result of this paper. The ratio Vc/c for singly quantized vortices 
are found to be between 0.16 — 0.31 for dipolar condensate with Edd = 0.15 and 0.15 — 0.28 for non-dipolar condensate 
in the range 6 > d > 0.5. As expected, the ratio Vc/c increases for doubly quantized vortices. 

In this paper, we have studied single and multiply quantized vortex dipoles in an oblate dipolar Bose Einstein 
condensate. We have shown that Vc is decreased while Vc/c is increased with increasing Edd in an oblate trap. The 
dynamics of vortex dipole in a dipolar BEC is worth studying. Helpful discussions with A. Kilic are gratefully 
acknowledged. 
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